Nematic Ferromagnetism on the Lieb lattice 
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We discuss the properties of possible ferromagnetic orders on the Lieb lattice. We show that 
the presence of a quadratic-flat band crossing point (QFBCP) at half filling will dramatically affect 
the magnetic ordering. In the presence of a weak on-site repulsive interaction, we find the ground 
state is a nematic ferromagnetic order with simultaneously broken of time-reversal and rotational 
symmetries. When the interaction strength increases, the rotational symmetry will restore at some 
critical value, and the system enters a conventional ferromagnetic regime. We also point out that the 
spin gap in both the nematic and conventional ferromagnetic phases is of the order of interaction. 
This observation suggests that these magnetic orders can be realized and detected in cold atomic 
systems with present technology. 
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The investigation on ferromagnetism is one of the cen- 
tral topics in condensed matter physics, and has at- 
tracted attention for nearly a century since the early age 
of quantum theory. Heisenberg showed that a system 
of localized spins would favor a fully polarized state by 
gaining exchange energy. |l| However, the same argument 
can not be simply applied in itinerant fermions, where the 
kinetic energy of the underlying system has to be consid- 
ered on an equal footing with the interaction effect. 

The discussion on the stability of itinerant ferromag- 
netism in a lattice system dates back to 1960's. Thou- 
less |2[ and Nagaoka |3j pointed out that a ferromag- 
netic ground state can be stabilized in any finite bi- 
partite lattice with infinite on-site repulsive interaction. 
Lieb showed that the stable region of ferromagnetism 
can be extended to arbitrary repulsive interaction, pro- 
vided that the number of sites are different for the com- 
posite sublattices. [4( The key ingredient in Lieb's argu- 
ment is the existence of a non-dispersive, or briefly fiat, 
band. When the flat band is partially filled, fermions 
tend to be spin polarized to minimize the interaction en- 
ergy, without any cost in the kinetic energy. In other 
words, since states in the flat band consists of localized 
Wannier functions, the ferromagnetic states can benefit 
from the exchange interaction as pointed out by Heisen- 
berg. Subsequent studies confirm the stability of ferro- 
magnetism in various models [I, Q, and generalize the 
idea to nearly- flat-band cases [7]. Experimental realiza- 
tion of the (nearly-) flat-band ferromagnetism has been 
proposed in a class of physical systems including atomic 
quantum wires |8|, q uantum-dot super-lattices [9], and 
organic polymers jlQj . 

Rapid progress in cold atom experiments has paved 
a new route toward the exploration of ferromagnetism 
in itinerant fermions. Thank for the extraordinary con- 
trollability of lattice potentials and interaction, several 
proposals has been made to realize (nearly) flat-band fer- 
romagnetism pj]-[l3|. In particular, Noda et al. inves- 
tigated two-component cold fermions loaded into a two- 
dimensional (2D) generalized Lieb lattice, and suggested 
that a ferromagnetic order can be stabilized in a wide 
parameter region. [l2| The 2D generalized Lieb lattice 
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FIG. 1: (Color online) (a) 2D Lieb lattice consisting of two 
square sublattices A and B. There are three sites (1, 2, 3) 
within a unit cell (dotted square), (b) The Lieb lattice can 
be realized by arranging three square optical lattices [lH ]: 
Vi = 2V [sm 2 (x/a) + sm 2 (y/a)}, V 2 = 2V [s'm 2 (x/2a + tt) + 
sin 2 (2//2a + 7r)], and Vj$ = Vb[cos 2 (x/2a + ?//2a) + cos 2 (x/2a — 
y/2a)]. Here, a is the lattice spacing, (c-d) Band structure 
of the Lieb lattice with Vb/t — 1. Note that the flat band 
degenerates with the upper dispersive band at the M-point, 
leading to a quadratic- flat band crossing point. 



consists of two square lattices (sublattice A and B), and 
has three sites per unit cell, as shown in Fig. [TJ With 
only nearest-neighbor hopping, the lattice topology sup- 
ports a flat band in the middle of two dispersive bands, 
hence can stabilize ferromagnetism when it is partially 
filled. 

Apparently, a sole fiat band does not exist in any real- 
istic physical system. In all proposed model lattices for 
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fiat-band ferromagnetism, the flat band is always asso- 
ciated with dispersive bands, and the effect of their ac- 
company is not fully understood. As a typical example, 
if the two composite sublattices of the Lieb lattice have 
the same depth, the flat band degenerates with the two 
linearly dispersive bands at the M = (ir, n) point. On the 
other hand, in general cases where the depth of sublattice 
B is shifted from that of sublattice A by an amount of 
T4, the flat band touches only one of the two dispersive 
bands at the M-point, as illustrated in Fig. [TJd), leading 
to a quadratic-flat band crossing point (QFBCP). 

In this manuscript, we show that the ferromagnetic or- 
der is dramatically affected by the existence of the QF- 
BCP. In the non-interacting level, the QFBCP is pro- 
tected by the time reversal (TR) and C4 rotational sym- 
metries. When a repulsive on-site interaction is present, 
since the 2D density of state (DOS) is singular, the QF- 
BCP becomes marginally unstable, leading to a sponta- 
neous broken of the TR and/or C4 rotational symmetries. 
We find that in the weak coupling limit the ground state 
is in a nematic ferromagnetic (NFM) order, character- 
ized by a broken of C4 point group down to C 2 . The 
C4 rotational symmetry will restore with increasing in- 
teraction, via a second or first order phase transition, 
depending on the value of Vb. Within a mean- field calcu- 
lation, we further map out the phase diagram at half fill- 
ing, and identify three phases including: (i) a semimetal 
with NFM order; (ii) a band insulator with NFM order, 
and (iii) a band insulator with conventional FM order. 
These magnetic orders have the potential to be realized 
and detected in cold fermions loaded in optical lattices 
with present technology. 

We consider spin 1/2 fermions loaded in the 2D Lieb 
lattice 

H = -t 4a C ^ + H^> i + E/^> W , (1) 

<i,j>,a i£B i 

where t is the hopping matrix, c\ a (ci a ) is the creation 
(annihilation) fermionic operator, Vb is the chemical po- 
tential offset for the sublattice B (i.e., the relative shift of 
the two sublattices), and U is the on-site interaction. In 
the context of cold atoms, V5 can be easily controlled by 
varying the relative intensity of optical lattices, and a re- 
pulsive interaction U > can be achieved and tuned via 
an adiabatic ramping to the upp er branch on the BEC 
side of a Feshbach resonance [15[. 

In the non- interacting case, the Hamiltonian can be di- 
agonalized in momentum space, leading to a band struc- 
ture consisting of three bands as shown in Fig.QJd). One 
of the three bands is completely flat, as required by the 
bipartiteness of the Lieb lattice. The flat band has Bloch 
wavefunction oc [0, — cos(fc £C ), cos(k y )] on the three sites 
within a unit cell, indicating the presence of local Wan- 
nier functions residing on the sublattice B and having 
opposite amplitudes between sites 2 and 3. When the 
chemical potential offset 14 = 0, the three bands are de- 
generate at the M = (n, tt) point, where the two linearly 
dispersive bands intersect with the fiat band. In general 



cases of Vb ^ 0, depending on the sign of V5, one of the 
two massive bands breaks the M-point degeneracy, and 
the flat band only touches the other dispersive band as 
illustrated in Fig. QJd). As a consequence, the M-point 
becomes a quadratic- flat band crossing point (QFBCP). 
Around this point, the effective two-band Hamiltonian 
reads 

where 5^ = k — M and 5k± = Sk x ± iSk y . From this 
effective Hamiltonian, it is clear that the cases of positive 
and negative chemical potential offset Vb are equivalent 
via a particle-hole transformation. Thus, we focus on 
systems with Vb > without loss of generality in the 
following discussion. 

The presence of a QFBCP at the M-point is the 
central feature of the Lieb lattice. This band cross- 
ing point (BCP) is protected by the time reversal (TR) 
and C4 rotational symmetries in the non-interacting case, 
and is characterized with a nontrivial topological index 
2tt. Such a putative topologically stable BCP becomes 
marginally unstable against infinitesimal repulsive inter- 
action U [16|, leading to a spontaneous broken of TR 
and/or C4 rotational symmetries, which drives the sys- 
tem toward a magnetic and/or nematic phase. Besides, 
since the BCP consists of a non-dispersive band, the infi- 
nite density of state (DOS) allows the possibility of filling 
the high momentum states with one single spin species 
without gaining any kinetic energy. As a consequence, if 
the on-site interaction U is repulsive, the system could 
easily favor a ferromagnetic phase for filling factors be- 
tween 1/3 and 2/3. 

Next, we focus on the case of repulsive Hubbard U > 
with half filling, and investigate the magnetic and ne- 
matic magnetic order within a mean-field (MF) level. 
The magnetic order is characterized by on-site magne- 
tization Bi and B2 for sublattices A and B, respectively. 
By minimizing the MF energy functional, we find that 
the system is stabilized by a staggered ferromagnetic or- 
der with Bi and B2 both oriented along the out-of-plane 
z-axis. In the weakly interacting limit [/Ct, the magne- 
tization on sublattice A is quadratically dependent on U 
with B lz w a 2 U 2 /9t, where a 2 = (2/tt 2 ) / d 2 k[cos 2 (fe E ) + 
cos 2 (k y )]/[V 2 /t 2 + 4cos 2 (^) + 4cos 2 (^)] 3 / 2 , and the 
magnetization on sublattice B is linearly dependent on 
the interaction with B 2z ~ — U/6. In the strongly cou- 
pling limit U ^> t, the interaction turns to be the only 
relevant energy scale, and both B\ z and B 2z are linearly 
dependent on U with B\ z ^ 2z = ±C//3, respectively. Here, 
we set Bi to be along the positive z-axis without loss of 
generality. 

To study the possibility of nematic magnetic order with 
both TR and C4 rotational symmetries broken, we allow 
the magnetization on sublattice B can be different for 
sites 2 and 3, and map out the zero-temperature phase 
diagram as shown in Fig. Efa). We find that the op- 
timized magnetization #$=1,2,3 are always along the z- 
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FIG. 2: (a) Zero temperature phase diagram for spin 1/2 
Fermi system on the Lieb lattice at half filling. The nematic 
ferromagnetic (NFM) order is favored for weak interaction, 
and is separated from the ferromagnetic (FM) order by a sec- 
ond order (thin solid) or first order (thick solid) transition line. 
Within the NFM regime, the system is a semi-metal (SM) for 
small U and V&, and undergoes a second order (thin dashed) 
or first order (thick solid) transition to become a band insu- 
lator (BI) with increasing U and V&. An example of the band 
structure for different phases (red dots) are shown in Fig. [3] 
(a-c). For two typical values of Vb/t = 1.5 and 2.3 (dotted 
lines), the magnetization on the three sites are shown in (b) 
and (c), respectively. 



axis, hence we consider only axial magnetic order in the 
following discussion. 

In the weakly interacting limit, the nematic FM phase 
is always favorable with an exponentially small magneti- 
zation difference \S\ = \B2 Z — Bs z \ ~ exp(— jt/U), where 
7 is a positive parameter depending on V& . This observa- 
tion is also confirmed by a perturbative analysis by treat- 
ing the magnetization difference as a perturbation to the 
existing FM order. When U is increasing from the weakly 
interacting limit, the NFM phase remains stable up to a 
critical value Z7 C , above which the C4 rotational symme- 
try restores and the system enters the conventional FM 
regime. For small values of V& < V c ~ 1.7, the phase 
transition between NFM and FM is of the second order, 
as identified by the condition V& = B\ Z — B2 Z . By increas- 
ing Vh > Vc, the NFM-FM phase boundary becomes of 
the first order, resulting from the competition between 
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FIG. 3: Band structure for (a) semimetal with NFM order, 
(b) band insulator with NFM order, and (c) band insulator 
with FM order. The solid (dashed) curves represent major- 
ity (minority) spin. These three cases are represented in the 
phase diagram Fig. [2fa) with red dots, and correspond to 
Vb/t — 2.3, U/t — 1.5, 5.0 and 8.0, respectively. The varia- 
tion of spin gap A s and band gap Ab for Vb/t = 2.3 is shown 
in (d). Note that the spin gap is always in the order of the 
interaction U. 



the corresponding metastable states. 

Within the NFM phase, the C4 rotational symmetry is 
spontaneous broken down to C2 by splitting the QFBCP 
into two Dirac points located along the direction of one of 
the principal axes for weak interactions. The two Dirac 
points have the same Berry flux 7r, in clear contrast to the 
case of graphene where the two Dirac points have Berry 
fluxes 7r and —it. In this case, the system is an anisotropic 
semimetal (SM) at half filling, with the Fermi surface 
shrinks to the two Dirac points as shown in Fig. [3fa). 
By increasing J7, the two Dirac points move toward the 
boundary of the Brillouin zone, and eventually disappear 
when the magnetization difference 5 — \B<i z — B^ z \ be- 
tween sites 2 and 3 exceeds the bandwidth of the first 
excited dispersive band. As a result, a full gap is open 
and the system becomes an anisotropic band insulator 
(BI) as depicted in Fig. [3fb). When the system enters 
the FM regime, the magnetization is large enough such 
that a finite gap is always present, and the system is a 
band insulator as shown in Fig. [3fc). 

Note that in both the NFM and FM phases, the spin 
gap remains in the order of the interaction strength U 
[See Fig. |3{d)]. This observation suggests that the mag- 
netic order on the Lieb lattice is robust against ther- 
mal fluctuations. Although the Mermin-Wigner theo- 
rem excludes the possibility of any 2D ferromagnetic or- 
der at finite temperature in the thermodynamic limit, 
the existence of such an order in a finite size system 
is perfectly allowed, provided that the coherence length 
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is comparable or exceeding the system size. Specif- 
ically, the coherence length of the magnetic order is 
£ ~ hcex.p(p s /T)/(kBT), where c is the spin- wave ve- 
locity at zero temperature, and the phase stiffness p s 
is in the same order of the MF transition temperature 
Tmf- Since the spin gap A s increases linearly with U, 
and the interaction can be tuned to be fairly large via 
an s-wave Feshbach resonance, the temperature required 
to observe the FM and NFM phases could be reachable 
within present technique. [17j 

The detection of the ferromagnetic orders can be im- 
plemented via an in-situ measurement [18|, LL9| , which is 
able to extract single site density distribution for differ- 
ent spin species, and hence the local magnetization Bi oc 
ni\ — n ii- If the system is prepared with equally pop- 
ulated two-component Fermi gas, we expect to see FM 
or NFM domains with opposite magnetizations, which 
are both resolvable within present experimental technol- 
ogy. Another possible detection scheme is to measure the 
single-particle dispersion with Bragg spectroscopy [20l or 
angle-resolved photoemission spectroscopy (ARPES) [21[ 
to extract the spin gap. 

In conclusion, we discuss the effect of a quadratic- 
fiat band crossing point (QFBCP) on the ferromagnetic 
(FM) order. Taking the 2D Lieb lattice as an exam- 
ple, we show that the QFBCP is marginally unstable 



against infinitesimal repulsive interaction, given an infi- 
nite density of state of the flat band. In the weak inter- 
acting limit, the ground state is a nematic ferromagnetic 
(NFM) order with time-reversal and rotational symme- 
tries broken. Within the NFM regime, the spontaneous 
generated magnetizations are different on sublattice B, 
and the QFBCP is broken into two Dirac points along 
one of the principal axes. In the strong coupling limit, 
the interaction U becomes the only relevant energy scale, 
and a conventional ferromagnetic (FM) phase with rota- 
tional symmetry restored is favored. We then map out 
the zero-temperature phase diagram within a mean-field 
analysis, and identify three regions including a semimetal 
with NFM order, a band insulator with NFM order, and 
a band insulator with FM order. We point out that the 
spin gap for all three phases is in the same order of in- 
teraction strength, which can be tuned via a Feshbach 
resonance. Thus, we expect these magnetic phases can 
be realized in two-component Fermi gases loaded in op- 
tical lattices at experimentally reachable temperatures, 
and can be distinguished via a species selective in-situ 
measurement. 
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